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Abstract 

We present the results of a global analysis of a class of models with an extended electroweak 
gauge group of the form SU(2) x SU(2) x U(l), often denoted as G(221) models, which include 
as examples the left-right, the lepto-phobic, the hadro-phobic, the fermio-phobic, the un-unified, 
and the non-universal models. Using an effective Lagrangian approach, we compute the shifts to 
the coefficients in the electroweak Lagrangian due to the new heavy gauge bosons, and obtain the 
lower bounds on the masses of the Z' and W bosons. The analysis of the electroweak parameter 
bounds reveals a consistent pattern of several key observables that are especially sensitive to the 
effects of new physics and thus dominate the overall shape of the respective parameter contours. 
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INTRODUCTION 



Despite the tremendous success of the Standard Model, there are still open questions 
that are unanswered and motivate further model-building. One of the most common model- 
building tools is to extend the gauge structure of the Standard Model. The simplest extension 
involves an additional U(l)x gauge symmetry (and thus an extra gauge boson Z'). One of 
the next-simplest extensions involves an additional SU(2), with the left-right model [Tj|2j[3j 
being perhaps the most widely-studied case of such models. On the other hand, given 
the extended gauge group SU{2)\ x SU (2)2 X U(l)x in the electroweak sector, there are 
many other models besides the left-right model that can be constructed, and these models, 
despite having a common fundamental gauge group, may have very different low-energy 
phenomenology. In this paper we present a unified, systematic study of many such models, 
which are commonly called (7(221) models in the literature. 

The most important feature of (7(221) models is the existence of new heavy gauge bosons, 
W and Z' . The existence of the gauge boson Z' has influences on the low-energy neutral- 
current processes, the Z-pole data at LEP-I and high energy LEP-II data [4j[5]. The exis- 
tence of the W boson has implications to the search of new physics beyond the Standard 
Model (SM) via studying charged-current processes. In low energy experiments, the most 
sensitive probes of charged currents come from flavor physics, such as the KK, bb mixing 
processes and semileptonic decays of the b quark [6] [7]. However, the low energy impact 
depends sensitively on the details of the flavor sectors, for which there is little experimental 
input [8] . There is thus a large uncertainty on the constraints on W' and its interactions. 

In this paper, we classify the G(221) models by the patterns of symmetry breaking sum- 
marized in Table [XT] (see section [Tl]) . Our main goals are to obtain the bounds on the masses 
of the W and Z' bosons for these various models, and, through the results of the global-fit 
analysis, to identify the key observables that are most sensitive to the new physics in these 
models. Our key results are that, at the 95% confidence level, the lower bounds on the masses 
of new heavy gauge bosons can be very light for breaking pattern I, which includes left- 
right, lepto-phobic, hadro-phobic and fermio-phobic models, for example, Mz> ~ 1.6 TeV 
and Mw ~ 0.3 TeV in the left-right model and hadro-phobic model; M z > ~ 1.7 TeV 
and My/i ~ 0.7 TeV in the lepto-phobic and fermio-phobic models. In breaking pattern 
II, which includes un-unified and non-universial models, because of the degeneracy of the 
masses of the W and Z', the lower bounds on their masses are quite heavy, for example, 
Mz' = My/i ~ 2.5 TeV in the un-unified model. 

We organize this paper as follows. In Section [ill we lay out the various G(221) mod- 



els and discuss the results of the relevant literature. In Section III, we give the effective 



Lagrangians, both at the electroweak scale (obtained by integrating out W and Z') and 
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below the electroweak scale (by integrating out the W and Z). In Section IV, we discuss 
the global-fit procedure and present our results obtained using the code Global Analysis for 
Particle Properties (GAPP) [9J, a software that utilizes the CERN library MINUIT |l(l and 
was used for the Particle Data Group global analysis [2E]- We also discuss which observ- 
ables are especially sensitive to the new physics contributions in these various models. We 
conclude in Section VI with a summary and outlook of our key findings. The Appendix 
contains the explicit effective Lagrangians for the G(221) models. 



II. THE G(221) MODELS 

We focus on the so-called (7(221) models having a S77(2)i x SU(2) 2 X U{l)x gauge struc- 
ture that ultimately breaks to f/(l) om . Relative to the Standard Model, these models have 
three additional massive gauge bosons, and their phenomenology depends on the specific 
patterns of symmetry breaking as well as the charge assignments of the SM fermions. For 
our studies, we consider the following different G(221) models: left-right (LR) [TJ|2][3], lepto- 
phobic (LP), hadro-phobic (HP), fermio-phobic (FP) P3JP2JP3], un-unified (UU) [H][l5], 
and non-universal (NU) [HI] [IZ] [IB] . The charge assignments of the SM fermions in these 
models are given in Table [T| and these models can be categorized by two patterns of sym- 
metry breaking (summarized in Table In]): 



• Breaking pattern I (the LR, LP, HP, and FP models): 

We identify SU(2)i as SU(2)l of the SM. The first stage of symmetry breaking then is 
SU{2) 2 x U{l)x — > U(l)y, giving rise to three heavy gauge bosons and Z' at the 
TeV-scale. The second stage is SU{2) L x U(l)y — > U(l) em at the electroweak scale. 

• Breaking pattern II (the UU and NU models): 

We identify U(l)x as U(l)y of the SM. The first stage of symmetry breaking is 
SU(2) 1 x SU{2) 2 SU{2) L . The second stage is SU(2) L x U(1) Y -> U(l) em at 
the electroweak scale. 

In addition to specifying the gauge group and the fermion charge assignments, a complete 
G(221) model should also include the ingredients of the Higgs sectors and the Yukawa 
couplings. While the observed relationships between the masses of W and Z bosons leave 
little freedom in the Higgs representation used for electroweak symmetry breaking (EWSB), 
we have freedoms in the choices of the Higgs representation used to break the fundamental 
G(221) gauge group to the SM gauge group. In breaking pattern I we assume the two 
simplest cases of symmetry breaking: via a doublet or a triplet Higgs. In the breaking pattern 
II we assume the simplest case of using a bi-doublet Higgs to achieve this symmetry breaking. 
The model-specific Higgs representations and vacuum expectation values (VEV's) are given 



in Table III For heavy Higgs boson, Wang et al. [19] used a non-linear effective theory 



approach to obtain an electroweak chiral Lagrangian for W . In our paper, by assuming a 
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TABLE I: The charge assignments of the SM fermions under the G(221) gauge groups. Unless 
otherwise specified, the charge assignments apply to all three generations. 



Model 


SU(2) 1 


517(2)2 


I7(l)x 


Left-rieht (LR) 




/ U R \ ( V R \ 


g for quarks, 
— 2 for leptons. 


Lento-nhobic (LP^ 






g for quarks, 
Ysm for leptons. 


Hadro-phobic (HP) 


M . ( ,/L ) 




Y^m for quarks, 
— \ for leptons. 


Fermio-phobic (FP) 


(;:)■(::) 




Kj M for all fermions. 


Un-unified (UU) 


(::) 


t) 


Y SM for all fermions. 


Non-universal (NU) 


fu L \ U\ 

V^i/ pt^nd \ eL J pt^nd 


fu L \ u\ 

V^/grd \ e i/ 3 rd 


Y SM for all fermions. 



TABLE II: Summary of the two different breaking patterns and the two different stages of sym- 
metry breaking in G(221) models. 



Pattern 


Starting Point 


First stage breaking 


Second stage breaking 


I 


Identify SU(2) 1 as SU{2) L 


SU{2) 2 x U{\) x -> I7(l)y 


Stf(2) L x C/(l)y -> U{l) em 


II 


Identify U{l) x as U(1) Y 


5C/(2)i x S£/(2) 2 -> 5C/(2) L 


SU(2) L X tf(l)y -> U{l) em 



light Higgs, we analyze the low-energy constraints by using a linearlized effective Lagrangian 
approach. 

The lepto-phobic, hadro-phobic, and un-unified models are, with the current set-up, 
incomplete because of gauge anomalies. It is entirely possible that the additional matter 
content used to address the anomalies can alter the low-energy phenomenologies and the 
results of our studies. Nevertheless, for completeness, we include these models with the 
current set-up in our studies, in which we focus on effects originated from the interactions 
of W and Z' bosons to the SM fields. In the cases of the lepto-phobic and hadro-phobic 
models, one can view them as transitions between the left-right models (where both right- 
handed leptons and quarks are charged under SU(2) 2 ) and the fermio-phobic model (where 
neither are charged). 

There have already been many theoretical and phenomenological studies of various (2(221) 
models, and we focus our brief literature review here mainly to those works that perform a 
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First stage 


breaking 










Rep. 


Multiplet and VEV 


LK-U, Lr-L) 
HP-D, FP-D 




(1,2,|) 








\ 






LR-T, LP-T 
HP-T, FP-T 




(1,3,1) 


$ = 


l 




1 

V2 


(- 


1 

n / 


UU, NU 




(2,2,0) 


$ = 


i 


v/2ir- 0° - irV 


1 

V2 


(: 


°) 

uj 








Second stag 


e breaking 












Rep. 


Multiplet and VEV 






LR-D, LP-D 
HP-D, FP-D 


£T~(2, 


2,0) 


H 












LR-T, LP-T 
HP-T, FP-T 


#~ (2,2,0) 


H 






V 






UU, NU 


#~ (1,2,1) 






C) 







TABLE III: These tables display the model-specific Higgs representations and VEVs that achieve 
the symmetry breaking of G(221) models. 



global fitting in the same spirit as our work. In the symmetric left-right model (where the 
couplings of the W are of the same strength as those of the W), Polak and Zralek obtained 
the constraints on parameters from the Z-pole data [20] and low energy data [2T] , separately. 
While for the non-symmetric case, Chay, Lee and Nam [22] considered phenomenological 
constraints on three parameters: the mass of the Z', the mixing angles (the analog of 
the Weinberg angle in the breaking of SU(2)r x U(1)x —> U(1)y) and the Z-Z' mixing 
angle £, by combining the precision electroweak data from LEP I (through e 1; 62, 63) and the 
low-energy neutral-current experimental data. For the non-symmetric case, the combined 
bounds at the 95% confidence level are 0.0028 < £ < 0.0065 and M z > > 400 GeV for all 0, 
while for the symmetric case, a more severe bound Mz> > 1.6 TeV is obtained. 

In the fermio-phobic model, Donini et al. [23J used the Z-pole and low-energy data, and 
the flavor physics data from flavor-changing neutral-current (FCNC) processes and b — > 57, 
to put constraints on the parameter space (W-W mixing angle a±, and Z-Z' mixing angle 
ccq) by fixing several sets of representative values of Mw' and x (strength of the coupling of 
the fermiophobic gauge group, relative to SU{2) L of the Standard Model). For the input 
parameters in the the ranges 100 GeV < M w > < 1000 GeV and 0.6 < x < 15, and for a 
low Higgs mass of 100 GeV, the best-fit values of | ckq | and |a±| increases with increasing x, 



6 



when holding Mw' fixed. On the other hand, when holding x fixed, increasing My/ 1 leads 
to an increase in the best-fit values of |«o| and a decrease in the best-fit values of |a±|. In 
the entire range of parameter space, the magnitude of the best-fit values of a and a± are 
at the percent level. 

In the non-unified model, Malkawi and Yuan [16] performed a global fit of the parameter 
space (x, 4>) using the Z-pole data, and found that the lower bound is Mz> = My/' > 1.3 TeV 
if no flavor physics is included. Chivukula et. al [21] used the data from precision electroweak 
measurements to put stringent bounds on the un-unified Standard Model [H] [IS] . They 
found a lower bound on the masses of the heavy W and Z' of approximately 2 TeV at the 
95% confidence level. 

III. THE EFFECTIVE LAGRANGIAN APPROACH 

To analyze the low-energy constraints, we will take an effective Lagrangian approach, 
and follow the general procedures laid out by Burgess in Ref. [25] to extract the effects of 
new physics. Although the details of each of these models are different, we first perform a 
generic analysis that can be applied to any G(221) model we consider in this work. 

Per the convention in Burgess [25], we denote the gauge couplings as g±, g 2 , and g x 
respectively for the gauge groups SU(2)\, S77(2) 2 , and U(l)x- The tilde (~) on the couplings 
and VEVs emphasizes the fact that these are model parameters, as opposed to quantities 
that can be directly measured in experiments, such as the physical mass of the Z boson. As 
an extension to the convention in Burgess |25j, we also denote with tilde (~) any combination 
constructed from the model parameters. We also abbreviate the trigonometric functions 



c x = cos (a;), s x = sin(x), and t x = tan(x). 



(1) 



A. Mixing Angles and Gauge Couplings 



We define the mixing angle at the first breaking stage as 




g x /g 2 (LR, LP, HP, FP models) 
h/h (UU, NU models), 



(2) 



and define the couplings 




(LR, LP, HP, FP models) 



(UU, NU models) 



(LR, LP, HP, FP models) 



(UU, NU models). 



(3) 
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The couplings g L and g Y are respectively the gauge couplings of the unbroken 577(2) L x U{1)y 
gauge groups after the first stage of symmetry breaking. Similarly to the Standard Model, 
for both breaking patterns we define the weak mixing angle (9) as 

t 9 (=twe) = ^. (4) 

For both breaking patterns, the electric charge (e) is given by 

1111 

e 2 91 gi 9 2 x 

and we also define a e = e 2 /A-K. 

With the angles 9 and 0, the gauge couplings can be expressed as 

„ je/(s§), (LR, LP, HP, FP models) 

9\ — < (6) 
1 e/(s^), (UU, NU models) 

- fe/(c^), (LR, LP, HP, FP models) 

#2 = < (7) 
1 e/(s^), (UU, NU models) 

- fe/(c S cj), (LR, LP, HP, FP models) 

= < _„ / (8) 



B. The Effective Lagrangian 

i. Gauge Interactions of Fermions 



In this sub-section we parameterize the gauge interactions of the fermions that is appli- 
cable to all the (2(221) models under considerations here. We will obtain both the SM-like 
effective theory applicable at the electroweak scale as well as the four-fermion effective the- 
ory below the electroweak scale. We do this by first building up the fundamental Lagrangian 
in stages, and then successively integrating out the massive gauge bosons. The Z-pole data 
measured at the electroweak scale, and measurements of the four-fermion neutral-current in- 
teractions are some of the most precise measurements to-date, and provide stringent bounds 
on new physics models. 

As discussed earlier, we consider the symmetry breaking to take two stages: 

517(2)! x SU(2) 2 x U(l) x -> SU(2) L x U(1) Y -> C/(l) em . (9) 

We denote the gauge bosons of the (7(221) models as: 

SU(2) 1 :W^,W^, 
SU(2) 2 : , W^, 

U(l) x : Xp. (10) 
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After the first-stage breaking, the neutral gauge eigenstates mix as follows 



W 6 = 



s^2% + ^ (LR, LP, HP, FP models) 
(UU, NU models) 

(LR, LP, HP, FP models) 
+ (UU, NU models) 

c ^ W ln ~ s 4> X n ( LR ; Lp ; HP, FP models) 



r?l — ) <P ~<p'~V \ ' ' -~ ' ~ * ~~l 

' - s^Wi^ (UU, NU models) 



and for the charged gauge bosons, we have 

w± ^ . (LR, LP, HP, FP models) 

" ~ ' SfWg, + CiW&, (UU, NU models) 

w , ± ^ J W*, (LR, LP, HP, FP models) 

" ^CiW^-SfW^. (UU, NU models) 

After the first stage of symmetry breaking, there is still an unbroken SU(2)l x U(1)y, which 
may be identified as the Standard Model gauge group. The gauge bosons W ±,z and B 
are massless, and only Z' and W f± are massive, with TeV-scale masses. The Lagrangian 
representing the heavy gauge boson masses has the form 

£ sta gc -i = m^z^z'i* + M^W'+W'-^, (13) 



VI 



where M\, and M^, are given in Table 

Before discussing the second stage of symmetry breaking, it is convenient to define, sim- 
ilarly to the Standard Model, (which will turn out to be the photon) and (approxi- 
mately the physical Z-boson) in terms of the massless gauge bosons and 

A, = (i-w*. + ^w^ + ^x, 

\9i 92 ^ 9x 

r3 



S § W* + c- R B u 



Z„ = cgWi - s § B„ (14) 



At the electroweak scale, the second stage of symmetry breaking takes place, breaking 
SU (2) x £7(1) — » U (l) cm . The Higgs vacuum expectation value (VEV) at the second stage not 
only gives masses to Z and W ± , but also induces further mixing among the gauge bosons 
W ± , Z, W' and Z'. The masses of the gauge bosons depend not only on the breaking 
pattern, but also on the group representations of the Higgs bosons whose VEV's trigger the 
symmetry breaking. For simplicity, for breaking pattern I, we consider only models with 
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either a doublet or triplet under SU (2)2, and do not consider models with both doublets 
and triplets. Introducing additional Higgses and VEVs would modify the masses of the 
W and Z' For breaking pattern II, since the first stage of symmetry breaking breaks 
SU {2)i x SU (2)2 to the diagonal subgroup, the masses of W and Z' are degenerate at this 
stage, and we only consider the case of an SU(2)i x SU (2)2 bi-doublet. For the convenience 
of typesetting, we also denote, for example, a left-right model with first-stage symmetry 
breaking triggered by an S77(2)-doub let (-triplet) as LR-D (LR-T). 

Although different breaking patterns and different group representations of the Higgs 
bosons will lead to different Lagrangians, we can write down the Lagrangian involving the 
gauge boson masses and fermionic gauge interactions in a general form 

C iund = l -M 2 z Z^ + 1 -{M% + AMf,)^'" + 5M 2 Z Z'^ 
+ W'+K + » + W'-R-v + Z'K°» 

(j, f.1 (j, 

+ W+J + » + W-J-» + Z fl J^ + A^, (15) 

where we have denoted the currents that couple to the primed gauge bosons (W' and Z') as 
K® and K^, and the currents that couple to the SM gauge bosons as J M , J° and J^. The 
SM-like currents have the familiar forms 

r = eJ2Q f 7Yf, (16) 
/ 

Jl = VW+M E i^U^p.k - sjQ f lu) , (17) 
/ 

J+ = ^= {d Llfl P L u L + e Llfl P L u L ) , (18) 
with an implicit sum over the three generations of fermions. The neutral currents (K®) and 



charged currents (K^), for the various models are summarized in Tables IV and |VJ We 
note the following features: 

• The residual SU(2)l x U(1)y is broken to the U(l) em , and there are now mass terms 
for the Z and W bosons, denoted as M zw . These masses have the familiar form 

*% = \(ft+fi)&, (19) 
Ml = \gy, (20) 

where the couplings g L and g Y are defined as in Eq. ^ for the two different breaking 
patterns. 

• There are mass-mixing contributions 5M Z w that induce Z — Z' and W — W' mixing. 



They are dependent on the breaking pattern and are given in Table VI 
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wy^u 






e7 M e 


LR 


1 ~ T-l 


(-|c^ 2 - \sfg x )P K 

1 ~ n 

-qS^ x P l 


(l c 4>92 + k s $9x)P R 

+ 2 S ^9xPl 


(-§c^ 2 + \s^g x )P n 

, 1 ~ n 


LP 


(\ c f92 - \sfg x )P n 
~\ s 4>9xP\ J 


(-\ c f92 - \sfg x )P K 
~I s 4>9xPl 


2 s <f>9x^L 


s ~n flp _i_ p > 
b 4>9x\2^ ' ftRj 


HP 


-sfg x (±P h + lP n ) 


s^x(\P,-\Pk) 


\2 C 4>92 T 2 b 4 > yx) Jr Yi 

+ \sfg x P h 


V 2 C 4>92 T 2 b 4,yx) r B.. 
+\ s 4>9xPi J 


FP 


-s-a^(-Pr + ~Pt>) 
q>9x \ g ftL i 


-s-a v (- Pr - - P ) 

q)9 X \ Q L 3 ftRj 


-s-a v Pr 
2 b (p9x^L 


s-o v (-Pt + P„) 


UU 










NU 


1 ( ) ft 


"3 ( ^ ) ft 


i ( ^ ) p L 


4 f c * 5 ! ) ft 



TABLE IV: This table displays the couplings g(f,f,Z') of the current K ^ = f^g{fJ,Z')f. 
For the top four models (LR, LP, HP, and FP), tan</> = g x /g2- For the lower two models (UU and 
NU), tancj) = 52/51- For the NU model (last row), the top values denote the couplings to the first 
two generations of fermions, and the bottom values denote the couplings to the third generation. 









LR 






LP 







HP 





^ 2 hPn 


FP 








UU 






NU 


1 ( %9x \ 





TABLE V: This table displays the couplings g(tp, f , W' + ) of the current K +tx = Jp^g®, £, 
For the top four models (LR, LP, HP, and FP), tan</> = g x /g2- For the lower two models (UU and 
NU), tan0 = 52/51- For the NU model (last row), the top values denote the couplings to the first 
two generations of fermions, and the bottom values denote the couplings to the third generation. 



There are additional contributions to the masses of the Z' and W' after the second 
stage of symmetry breaking, which we denote as AM§/ W i. They are also dependent 



on the breaking pattern and are given in Table VI 
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M 2 , 


M 2 ,, 


AM 2 , 


AM 2 ,, 


5M 2 


5M 2 , 


LK-D, Lr-u 
HP-D, FP-D 


\(gl + ~9 2 x )ui 


1 nV, 2 


c 2 

4> 7.2-2 

-$-92 v 




r 2 

Ay ~ ~ ~ ~2 

-j-gmgxv' 


-\gig2v 2 s 2 p 


LR-T, LP-T 
HP-T, FP-T 


(g\ + 9 2 x)ul 


\g\^T 


4 

-fg2 v 




cj 

-^hMxv 2 


-\gig2v 2 s 2 - p 


UU, NU 


\(~g 2 + ~g 2 )u 2 


\{g\^g\)u 2 


4 

4> 727.2 

~^g2 v 


4 

^g2 v 


4 

-^gig2g x v z 


-\~gihv 2 sl 



TABLE VI: This table displays the model-dependent parameters M z , w , in Eq. (13), and AM Z , w , 



and 5M ZW in Eq. (15). 



Therefore, the gauge boson mass terms can be written as 
r 

' m ass 



w + W H 



M 2 W 



SM 2 , 



5m 2 , m*,, + am;, 

/ 

M| 6M Z 
\ 5M 2 M 2 , + AM 2 , J 



( A\ 



(21) 



In Table III, we expect that the scale u 2 of the first-stage breaking is much larger than the 
electroweak scale v 2 . We work to leading order in v 2 /u 2 , and so if we take the approximation 



M 2 , w , > M 



ZWi 



6M 2 ,AM 2 , 



(22) 



we can expand in large M|, w ,. To order 0(My^, z ,), the mass eigenstates, denoted without 
the hats ("), are given by (similarly for the charged gauge bosons): 



SMI a, 

2^' 



M 2 , - M z 



Z' 



5M 2 Z 



Mf , - Mf 



Z' 



(23) 
(24) 



Now we can rewrite the fundamental Lagrangian in terms of the mass eigenstates for both 
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neutral and charged gauge bosons 



/'•mass 
^fund 



Ml - -^A 



z„z» - [ i\4 - 2^ 



5M 



Mf, 



' 1 M z> 



+ A (t J» 



Ml a r 2 



M, 



+ (+<-► 



IF' 



(25) 



We can now obtain the effective Lagrangian by successively integrating out the massive 
gauge bosons. In the basis of the mass eigenstates, integrating out Z' and W (whose masses 
are expected to be at or above the TeV scale) results in an effective Lagrangian valid at the 
electroweak scale: 



cff 




<5M| 
M% 



Z^ + 



M 2 W 



8Mb 



w +fl w: 



+ z ^ | 



Ml 



K 0lM K c , 



2M% M 



M 2 ^ ' * 

1V1 W' 



(26) 



From Eq. (26), we see that the low-energy effects of the heavy gauge bosons are parame- 



terized by the shifts in the masses of the W and Z gauge bosons, and in the shifts of their 
couplings to the fermions, and additional four-fermion interactions. 

We can further integrate out the Z and W ± gauge bosons (again to leading order in 
M^> z 1 ) • We then have the four-fermion interactions 



cff' 



2M| 



jo** + Ml 

Ml, 



M 2 W 



M, 



w> 



+K+K-*)] 



r4 " M 2 » » 



Mi 



<)M*r J+J ^_6M^ 



M 2 W 



(27) 



Before we can compare the predictions of Eq. (27) with experimental results for the different 



(7(221) models, we first have to properly define some experimental input values (for example, 
the Fermi constant Gp) for the G(221) models under study. We will discuss this in Section 

EH 
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2. Triple Gauge Boson Couplings 



In the basis denned through Eqs. (11), (12) and (14), the triple gauge boson couplings 



(TGCs) g(ZW + W~) and g(AW + W~) have the standard forms 

g(ZW + W~) = -g L cos 9, (28) 
g(AW + W-) = e. (29) 

In the basis of mass eigenstates, however, we expect there to be a shift to these couplings 
because the mass eigenstate Z (W) now is a mixture of Z (W) and Z' (W). However, 
because of QED gauge invariance, the AW + W~ coupling does not receive a shift. On the 
other hand, the ZW + W~ coupling does shift, and we shall discuss in turn this shift for the 
two breaking patterns. 

In breaking pattern I (LR, LP, HP, and FP models), in the hat (") basis of the gauge 
bosons, the Lagrangian contains ZW'W' and Z'W'W' vertices in addition to the typical 
ZWW vertex. Since the overlap between W' and the light mass eigenstate W is of order 
0(Mw>), contributions from g (zW'W*) and g ( Z'W'W') to g (ZWW) are at least of order 

O(M^). As we are only working to leading order in O(M^), there is no shift due to these 
additional interactions at this order. 

For breaking pattern II, the story is similar. There are no ZWW 1 nor Z'WW vertices, 
only ZW'W' and Z'W'W interactions. The contributions to the ZWW coupling are sup- 
pressed by fourth powers of the heavy masses z ,, and thus of higher order than those 
kept in the effective theory. 

In both breaking patterns, however, there will be a shift to the ZWW-vertex due to a 



shift in 9 (cf. Eq. (47)) , the counterpart of the Standard Model weak mixing angle 9 , 
as defined in our fitting scheme. The LEP-II experiments, however, do not directly probe 
the ZWW-vertex, but instead infer the ZWW-vertex through the process e + e~ — > W + W~ 
assuming SM couplings for all other vertices. To properly compare the relationship between 
the experimental measurement of the ZWW-vertex and the theoretical shifts in the (7(221) 
models, we would have to take into account all the other shifts in the couplings that enter 
the process e + e~ — > W + W~. We will discuss this in further detail in Section fv} 



3. The Yukawa and Higgs Sectors 

We complete our discussion of the effective Lagrangians of the (7(221) models with a 
brief discussion of the Higgs sectors and the Yukawa interactions. It is important to stress, 
however, that despite the complexity of the Higgs sectors and Yukawa interactions, our 
results of the global analysis only depend on the gauge interactions of the fermions, and 
not on the details of the Yukawa interactions. This is because we work only with those 
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observables involving gauge interactions (which excludes, for example, the branching ratio 
Br(6 — > 57)), and keep only tree-level contributions originated from the new physics. 

We discuss the Higgs sectors of the two breaking patterns separately. In breaking pattern 
I, we take as an example the left-right model, where the electroweak symmetry is broken by 
a bi-doublet (LR-D). This is necessary because the VEV's of the bi-doublet should generate 
the fermion masses, and the right-handed fermions now are doublets under the SU (2)2. 



With the bi-doublet H in Table III, we may have Yukawa couplings (similarly for leptons) 
such as: 

-C D Q R (y Q H + y Q H) Q L + h.c, (30) 

with H = —ir2H*T2, and where 3^q and 3^q have flavor structures that may be related by 
imposing additional symmetries (for example, left-right parity) on the model. In any case, 
unlike the Standard Model where we can solve for Yukawa couplings in terms of fermion 
masses and the Higgs VEV, in G(221) models there are more free parameters in the Yukawa 
sectors. These parameters can lead to interesting flavor phenomena, particularly in the 
arena of neutrino physics, and have been studied in detail in the literature (see, for example, 
Mohapatra et. al. [27J). On the other hand, the details of the Yukawa sectors do not affect 
the gauge couplings of the fermions at leading order and therefore do not affect the results 
of our analysis. 

In breaking pattern II, in addition to those Higgs bosons that are required to break the 
electroweak symmetry, it may be the case that the Higgs sector needs to be extended to 
generate fermion masses. This is because, with the current set-up, the Higgs boson that 
generates EWSB can couple only to leptons (in the case of un-unified model) or fermions of 
the third generation (in the case of the non-universal model). With additional Higgs bosons, 
the structure of the Higgs potential may mimic that of the two-Higgs doublet models. Again, 
as with breaking pattern I, there are more degrees of freedom than can be determined from 
the fermion masses, but the details of the Yukawa interactions do not affect the results of 
our paper, which only depend on the fermionic gauge interactions. 



IV. THE GLOBAL FIT ANALYSIS 



In this section we illustrate our procedure for performing the global-fit analysis to obtain 



constraints on new physics contributions. From Tables III and IV, we see that the G(221) 



models contains six (five) parameters for the first (second) breaking pattern: three (two) 



VEV's {u D)T , v sin f3, v cos f3} in Table III and three gauge couplings {<7i, <72> fe} i n Table IV. 



(For breaking pattern II, there are only two VEV's {u, v}.) Compared to the gauge sector 
of the SM, which contains only three parameters (two gauge couplings and one VEV; g L , g Y 
and v), there are three (two) additional parameters, and our goal is to: 
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• find a useful parameterization of these three additional parameters so as to parame- 
terize the effects of new physics, and 

• determine the constraints on these parameters from electroweak precision measure- 
ments through a global-fit analysis. 

We discuss these two steps in detail in turn. 



A. Parameterization 



As stated above, the (7(221) models contain six (five) parameters in the gauge sector: 

{<?i, h, 9x, u D (u T , or u), v 2 , /?}, (31) 



where the parameter (3 only exists in models with breaking pattern I. Using Eqs. (|6j),([7j), 
and (j8J), an equivalent set of parameters is 

{a e , 6, 0, x, v 2 , s 2 ^}, (32) 

where x is defined as 

for LR-D, LP-D, HP-D, and FP-D 

for LR-T, LP-T, HP-T, and FP-T (33) 
for UU and NU. 

As we expect x to be large (x > 100), we work to leading order in x~ l . 

In addition to these parameters, the loop-level predictions will require the values of the 
masses of the top quark (mt) and the Higgs boson (M#). For each (7(221) model, we perform 
two separate analyses with regard to these parameters. In one analysis, we fit these two 
parameters, m t and M#, in addition to the model parameters. In a second analysis, we fix 
these two parameters at the best-fit SM values. 




With regard to the parameters in Eq. (32), we will take three reference observ- 
ables to constrain three combinations of the parameters and perform a global-fit over 
{x, (j), s 2 s, fn t , Mil). The bar ( - ) over m t indicates that we will use the top quark mass as 
defined in the MS-scheme. We take as reference observables the experimental measurements 
of 

• the mass of the Z boson (Mz = 91.1876 GeV), determined from the Z-line shape at 
LEP-I. 

• the Fermi constant (Gf = 1.16637 x 10~ 5 GeV -2 ), determined from the lifetime of the 
muon, 

• the fine structure constant (a" 1 = 127.918 at the scale M z ). 
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Our task then is to express the model parameters, cf Eq. (32) 



{a e , 9, v , x, 0, s 2/ 3, m t , M H }, 
in terms of the reference and fit parameters 

{a e , M z , G F , x, <j>, s 2j g, fn t , M H }. 
That is, we want the relationships 

model parameters reference parameters fit parameters 



(34) 



{a e , 9, v , x, 0, s^,m t ,M H } <^> { a e , M z , G F , 



x 



s 2 s, m t , M H } (35) 



Since {x,(fi, (3, m t ,M H } appear in both the model and fit parameters (by construction), 
we only have to solve for {a e ,8,v 2 } in terms of the reference and fit parameters. This can 
be done by analyzing how the reference parameters are related to the model parameters. 

1. Electric Charge 



The electric charge in the G(221) models is the gauge coupling of the unbroken i7(l) em 
group, which we have parameterized as e in Eq. Q. There are no tree-level modifications 
to the wavefunction renormalization of the photon, so we then simply have the relationship 



(36) 



2. The Fermi Constant 



The Fermi constant, Gf, is experimentally determined from the muon lifetime as [28] 

,2 



T, 



= G 2 F ml 

1927T 3 



ml 



, 771 „ 
1+0 ' 



M 2 W 



1 + 



1 



167T 1 



where the precise forms of the higher-order corrections are given in Ref. 
these higher-order corrections, the SM contribution to the muon lifetime is 



d m 5 

192-32tt 3 M^ 



and, using the SM relation 4M^/ = gjv , we obtain 



Gj 



1 



(37) 
Neglecting 

(38) 



(39) 
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In the G(221) models, we have extra contributions to the four-fermion charged-current 



effective theory below the electroweak scale, cf Eq. (27), 



rCCAS 



M 2 W 



-J + J- 



M 2 W , 



K + K~ 



M 2 W 



(J + K- + K-J + ) + 



and these contributions will modify the SM relation in Eq. (39). In principle, the fermionic 



contributions to K^ can have both left- and right-handed components and differ among 
the different generations. However, for the G(221) models we consider here, Ifr couples 
universally to the first two generations. Furthermore, K^ is either purely right-handed (the 
LR, HP, LP, FP models) or purely left-handed (the UU and NU models). We therefore focus 
on these special cases instead of performing the general analysis. 

We first consider the case that is purely right-handed. The contributions to the 
amplitude come from J J, JK, and KK operators that do not interfere with one another in 
the limit of neglecting the masses of electrons and neutrinos. The squared-amplitudes from 
the JK and KK operators are of order 0(M^,) ~ 0(x~ 2 ) at leading order, and we do not 
keep these contributions. The Fermi constant is then given by 



G p 

71 



~9l 



8M^ 



1 + 



(for breaking pattern I), 



(40) 



independent of the details of K^. The expression of Gp, which depends on the details of 
the Higgs representation, is written in terms of model parameters as 



Gi 




(for LR-D, LP-D, HP-D, and FP-D) 
(for LR-T, LP-T, HP-T, and FP-T) 



(41) 



Though the left-right and right-right current operators do not contribute to the total 
muon decay rate at the order C(x _1 ) , they do contribute at leading order to the Michel pa- 
rameters (for a detailed discussion of the Michel parameters, see the Muon Decay Parameters 
article in the Particle Data Group (PDG) 



In the case that is purely left-handed, all the charged-current operators in Eq. (27) 
contribute, and Gp is given by 



Gp 

71 



~9l 



8M^ 



1 + 



M 2 W 
M 2 WI 



g 2 w , JM&gw. | 6M& 



~9l 



M? 



w 9l 



(for UU and NU) 



(42) 



where gw' can be looked up in Table |V| For the UU and NU models, these contributions 
cancel each other, and we are simply left with 

1 



G f 



V2 



r.2 



(for UU and NU). 



(43) 
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We can rewrite our results in a more suggestive manner by defining the SM VEV (v 2 
without tilde ~) through the Fermi constant 



We then have 



v 2 1 + -S 



2d 



V 2 1 + 



X 

s 2 ~ 
2x 



\/2<V 



(44) 



(for LR-D, LP-D, HP-D, and FP-D) 



(for LR-T, LP-T, HP-T, and FP-T) 



(45) 



v 2 . (for UU and NU) 



3. Z-Mass 



In our effective theory approach, the mass eigenvalue of the Z-boson is given by (using 



Eq. (25), Table VI, and a e = a e ) 

5M% I general form from the 



Mi = Mi 



M 7t I fundamental (7(221) Lagrangian 



r,2 



a e 7TV 

s\c\ 



i 

X 



(for LR-D, LP-D, HP-D, and FP-D) 



■0- ( 1 - ^ , (for LR-T, LP-T, HP-T, and FP-T) 



a e Tcv 



(46) 



a e 7iv 2 



1 - -± ] , (for UU and NU) 



Solving Eq. (46J) for c~s~, and using Eqs. (|41|) and (|43|), we can solve for 9 in terms of the 
reference and fit parameters 



2 2 , 

SnC„ = < 



s 2 2 



„2 2 



„2 2 



1 £ (4 S 2/3 



-if 1 
X I 4 



C 4 - - V- 



(for LR-D, LP-D, HP-D, and FP-D) 
(for LR-T, LP-T, HP-T, and FP-T) 



1 - 



x 



(47) 



(for UU and NU), 



where 9 (without a tilde ) is defined in terms of the reference parameters 



sin 9 cos 9 = 



TXOLf, 



V2M 2 G F 



(48) 



Eqs. (36), (45), and (47) then enable us to translate all the model parameters to reference 



and fit parameters. 
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B. Corrections to Observables 



In this subsection we illustrate the corrections to several example observables that we 
include in our global analysis. These examples elucidate the procedures we had outlined 
earlier, and we will refer to these results when we discuss the observables included in our 
global analysis. 



1. The Z-Partial Widths T(Z ff) 



As a first example, we can then consider the Z — >■ // partial width, which at tree-level 
has the expression in the Standard Model 



r(z -> //) 



12tt 



Mz {g 2 v + 9a) , 



where n c = 3 if / is s quark, and n c — l for leptons, and 



9v = 



2s e c e 



T 3 / L -2Q'sin 2 #), 



9a 



2s e c e 3L ' 



(49) 

(50) 
(51) 



where T^ L and Qf are respectively the weak-isospin and electric charge of the fermion /. 



In the (7(221) models, the partial decay width can be written in terms of model parameters 



as 



r(z //) 



12n 



5M A Z 
2M 2 Z M 2 Z , 



[9 z v {f)f + [9 z Af)f) , (52) 



where 8M Z , M z ,,g^(f), and g^(f) depend on the details of the model. For models that 
follow the breaking pattern I (LR-D, LP-D, HP-D, FP-D), the couplings have the form (to 
order O^ 1 )) 



2snCn 



c 2 7 



2x 



T( R cl - (X[ + X> R )s 



FAf) = 



Tl R c\ - (X' R - Xl)s 



A„2 



(53) 
(54) 



where X[, and X R , and T 3 ^ are respectively the left- and right-handed fermion charges under 
the U(l)x, and the ^-component isospin under the SU(2) 2 (which is identified as SU(2)r in 
left-right models). Expressing 9 in terms of the reference and the model parameters through 
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Eq. (47) and collecting terms of 0(x 1 ), we have (in units of GeV) 



T(Z -> //) = T(Z -+ ff) m 

a\ (-0.446 (Q f ) 2 + 1.773 Q f T( L - 0.310 Q f T? 



x 



?,R 



3L.) 



-0.310 Q f X f R - 0.664 (Tj 
+ sf (0.582 (Q 7 ) 2 - 1.91 + 0.620 Q J T( R + 0.310 Q f X 

+ s\ B (0.136 (Q / ) 2 - 0.136 Q f Ti L - 0.664 (T^ 
-0.136 (Q 7 ) 2 + 0.136 Q f T[ L - 0.310 Q / T 3 / JJ + 0.664 (T? L ) 2 
(for LR-D, LP-D, HP-D, and FP-D) 



(55) 



where T(Z — > //)sm is given by Eq. (49), and we have used the numerical values of the 
reference parameters. 



2. The Mass of the W-boson 

As a second example, we compute the mass of the VT-boson in the G(221) models. The 
SM expression, for the same set of reference parameters {a, Mz, GV}, is given by 



M w = M z c e , 



(56) 



where 6 is defined in terms of the reference parameters in Eq. (48). In the G(221) models, 
the mass of the l^-boson has the general form 



M W = M W \1 



2M 2 ,M 2 / , 



(57) 



More specifically, in terms of the model parameters for the individual models, we have 



M w = { 



2F I 1 - 2§ J ( for LR - D ' Lp - D ' Hp - D ' Fp - D )' 

Sg ( 1 ~ if j ( fOT LR " T ' LP " T ' HP " T ' FP " T )' 

Iff l-g) (forUU,NU). 



(58) 
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Using Eqs. (36), (45), and (47), we can convert all the model parameters to reference and 



fit parameters 

M z cos 9 



Mi 



w 



M z cos 9 
M z cos 9 



1 + 



1 + 



1 + 



2x a - si \°i s w 



C — Sn 



c 4 s 2 
4> 20 



(for LR-D, LP-D, HP-D, FP-D), 
(for LR-T, LP-T, HP-T, FP-T), 



(for UU, NU). 



(59) 



C. Implementation of the Global Fit and List of Observables 

For a measured observable O cxp , the SM prediction can be broken down into the tree- 
and loop-level components 



th 



Of^{m u M 1 



n) 



(60) 



where th is expressed in terms of the reference parameters. Since the top quark mass (fn t ) 
and the mass of the Higgs boson (Mh) enter into the loop-calculations in the SM, a global 
analysis of precision data and direct detection data can be used to constrain Mh- In the 
G(221) models, we can express the theoretical prediction as 



= Of^ c + Of^(m t , M H ) + O^ix, 4>, £), 

where 0*p trcc is of the order (9(1/5;), and we assume that 

1 



x 1 ~ 



167T 2 



~ O 



th.loop 

SM 



(61) 



(62) 



That is, the Born-level new physics contributions from the G(221) models are numerically of 
one-loop order, and loop corrections involving new physics are of two-loop order O ( jq 1 2~ ), 
which we discard in our analysis. 

To compare with precision data (from LEP-1 and SLD) and low-energy observables, we 
calculate the shifts in observables 0*p trcc (x, 0, /?), as in the previous examples of the partial 
decay widths of the Z-boson and the mass of the PF-boson, and we adapt these corrections 
into a numerical package GAPP [9]. GAPP then computes Og^ trcc and Of^ oop (m t , M H ) jl2] , 
together with the 0^p tiee (£, <f>, /3) to find the best-fit values of the fit parameters and the 
confidence level contours using the CERN library MINUIT [TTJj . 

We perform a global fit over the following classes of observables 

• LEP-I Z-pole observables: the total Z- width (r^), left-right asymmetries (A LR ), and 
related observables. 



the mass (M w ) and decay width (T w ) of the PF-boson, 
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• the tau lifetime r T , 



• the ratios of neutral-to-charged current cross sections measured from neutrino-hadron 
deep-inelastic scattering (DIS) experiments (R u = cr^/a^ and similarly defined for 



• effective vector and axial- vector neutrino-electron couplings (#y e and g 1 ^), 

• weak charges (Qw) °f atoms and the electron measured from atomic parity experi- 



Detailed information on these observables can be found in PDG [2H], and here we only briefly 
summarize the observables. The set of the observables included in our analysis is the same 
as that used in the PDG analysis [28], with two exceptions. 

• First, we do not include the anomalous magnetic moment of the muon and the decay 
branching ratio b — > s'j. At leading order, these observables are of one-loop order, and 
they depend on the details of the extended flavor structure of the (2(221) models. In 
this work, we assume W bosons only couple to fermions in the same generation. 

• Second, we include the measurements of the decay width of the IV-boson, which are 
not included in the PDG analysis. However, because of the comparatively low precision 
of these measurements, this observable turns out to be insensitive to the new physics 
contributions from the (2(221) models. 

In total, we include a set of 37 experimental observables in our global-fit analysis. 

Before we give a brief discussion on each of these classes of observables, we note that for 
some low-energy observables, such as the measurements from the atomic parity violation 
and neutrino-neucleus DIS experiments, we implement the shifts in the coefficients of the 
relevant four-fermion interactions, and rely on GAPP to compute the theoretical predictions 
based on these modified coefficients. The expressions of the coefficients of the four-fermion 
interactions are given in the Appendix. 

For the ease of typesetting in the following subsections, we introduce the abbreviation 
for the various forms of the fermionic currents 




ments. 



(/ l/2 )^ = / l7 ^(l- 75 )/, 
(/i/aJ^/irCl+Ti)/: 

(hf 2 Y v = hrf2, 
(hf 2 y A = / l7 ^ 5 /2. 



2- 



2- 



(63) 



23 



1. Precision Measurements at the Z-Pole 



The precision measurements at the Z-pole (including LEP-1 and SLD experiments) fall 
into two broad classes: observables that can be constructed from the partial widths (for 
example, in Eq. (55)) and the asymmetry (constructed from the couplings in Eq. (50) and 
(51)). We discuss these two classes in turn. 

In addition to the total width Tz, there are also the following measurements: 

12n 

a - d = jpf2 ■ T z (e~e + ) T z (had.) , (64) 

z z 

m = T ^W' for£ = e,/i,r, (65) 

R{q) = J Z } q ^ \ , for q = u,d,c,s,b, (66) 
T z (had.) 

U{S) = R(u) + R(d) + R(sY (67) 
where Tz(ff) is the partial decay width T(Z — > ff), and 

T z (had.)= Tz ^- ( 68 ) 

q=u,d,c,s,b 

The left-right asymmetry Al.r(/) is defined as 

b? CO] 2 + \s%U)f 

where (/) and g^{f) are the couplings of the fermion / to the Z-boson: 

c d z^uV.rh + gi(f)fn^U). (70) 

From the quark branching ratios R(q) defined above, the hadronic left-right asymmetry Qlr 
can be defined as [H] 

Q LR = R (9)A LR (q) - R(<l) A LR{<i)- (71) 

q=d,s,b q=u,c 

A second class of asymmetries, the forward-backward asymmetries Afb{J), emerges from 
the convolution of the Al.r(/) asymmetries with the polarization asymmetry Ai R (e) of the 
electron. The hadronic charge asymmetry Qfb is defined accordingly [9] [29] 

A FB (f) = 3 -A LR (e)A LR (f) } (72) 
Qfb = ^A LR (e)Q LR . (73) 
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2. The Tau Lifetime 



In terms of model parameters, the expression of the tau (r) lifetime is similar to the muon 



(fi) lifetime in the (7(221) models, cf. Eq.(38), with the obvious replacement of in the 



fx lifetime by m T in the r lifetime. This is true even in the non-universal (NU) model, in 
which third generation fermions transform under a different gauge group compared to the 
first two generations. In the four-fermion effective theory of the NU model, only interactions 
involving two pairs of third-generation fermions receive new physics contributions, and the 
interactions involving one pair of third-generation fermions with one pair of light-flavor 
fermions (those responsible for the decay of the r) are the same as those between two pairs 
of first two generations of fermions (those responsible for the decay of ft). This is similar 
to the case of the un-unified model, where only interactions involving two pairs of quarks 
(qq)(qq) receive new physics contributions, while the (qq){££) interactions are the same as 
the {££){££). The lifetime r T can be calculated at tree level as 



_ x G F m\ ( ^ o m; 

Ml' 



[ 1 + 3 W >• (74) 



in the SM. The dominant new physics contribution from G(221) models can be captured in 



the shift of M w as shown in Eq. (59). 



3. vN Deep Inelastic Scattering 

The uN deep inelastic scattering experiments probe the coefficients e L (q) and e R (q) (for q 
being u or d) that parameterize the neutral current Vvqq interactions below the electroweak 
scale 

q=u,d 

The DIS experiments measure the ratios of neutral-to-charged current cross sections 

Rv = Vvn/VvNi R v = a vNl a vNi ( 76 ) 

which can be written in terms of e L (q) and e R (q) as 

R v = (1 - 5) [a L {u)e 2 L (u) + a L (d)e 2 L (d) + a R {u)s 2 R (u) + a R {u)e 2 R (d)] , (77) 
Rp=(1- 5) [a L {u)s 2 L (u) + a L (d)e 2 L (d) + a R {u)e 2 R (u) + a R {u)e 2 R (d)] . (78) 

The coefficients 5 and a L)R are related to the nuclei form factors that are experiment specific. 
These coefficients are included in GAPP, and we implement only the corrections to e L (q) 
and e R (q). 
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4- ve Scattering 



The most precise data on neutrino-electron scattering comes from the CHARM II [30] 
experiment at CERN that utilized and V^. The relevant parameters e L (e) and e R (e) are 
defined similarly as in the uN scattering 



C D 



G 



F /_ 



V2 



( e ) Hl+ £ - ( e ) N3- 



(79) 



We can further define 



2v 



(e) + e L (e) , 
£r (e) - e L (e) , 



(80) 
(81) 



which are related to the measured total cross sections cr^ e c and cr^ e c or their ratio 0"^ e °/c"^ e c . 
In the limit of large incident neutrino energies, E v 3> m e , the cross sections are given as 



a 



G\vn e E v 

G 2 F m e E u 
2^ 



{97 + 97? + \{97 ~ 97? 
(97 -97? + \(97 + 97? 



(82) 
(83) 



We implement corrections to the couplings due to new physics in GAPP and compute the 
cross sections that are used in the global-fit analysis. 



5. Parity Violation Experiments 

We consider observables from three different measurements: atomic parity violation 
(APV), M0ller scattering (e~e~ — > e~e~) [21], and eN DIS. These experiments measure 
the weak charge (Qw) of the electron [31], caesium-133 [32] [33] and thallium-205 nuclei 
[34j|35j. Before defining the weak charge, it is useful to parameterize the coefficients of the 
(ee)(qq) and (ee)(ee) interactions in terms of Ci q , C2 q , and C± e as 



J4) 



The weak charges of the quark and electron are defined as 

Q w (q) = 2C lq , Q w (e) = 2C le . 



(85) 



We can express the SM tree-level couplings of quarks to the Z-boson as C D Z^J^ , where 



J*=\9l(q)\- Qw(q)(qq) v +(qq) 



A,H 



(86) 
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and the ± on the axial- vector term is the opposite sign of the T^ q . Hence Qw{q) can be 
interpreted as the ratio of the vector current to axial- vector current coupling of quark q to 
the Z-boson: 

g z (q) 

Qw,sm{q) = | V T , \ I ■ (87) 

\9l(q)\ 

The weak charges of the nucleons and nuclei can be built up from those of the quarks 

Qw(p) = 2Qw(u) + Qw(d), (88) 
Qw(n) = Qw(u) + 2Q w (d), (89) 

and for nucleus A Z (with atomic number Z and mass number A), which contains Z protons 
and N(= A — Z) neutrons, 

Q w ( A Z) = Z ■ Q w {p) + N ■ Q w {n) (90) 
= 2 [{Z + A) ■ C lu + (2A - Z) ■ C ld \ . (91) 

There are also measurements of certain linear combinations of the coupling coefficients 
Ciu and C\d from polarized electron-hadron scattering data [SB]. The particular linear 
combinations, determined by the experimental data, 

C x = %C Xu + 4C xd , 

C 2 = -4C lu + 9C W , (92) 

are included in our global analysis. 



V. RESULTS 

A. Global Analysis 

In this section, we present the allowed regions of parameter space based on the global-fit 
analysis. A testament to the success of the SM is that, for all the (7(221) models, the global 
fitting pushes x to large values, decoupling the effects of the new physics. This is presented 
in Figs, [l] and [2j where we show the 95% confidence level (C.L.) contours on the x — 
plane. 

In addition to the constraints from the precision and low-energy data, we also require 
cos0 (sin0) to be greater than 0.1 (0.18) for the first (second) breaking pattern so that all 
the gauge couplings in Eq. Q and & are perturbative and do not exceed \/47r. These 
constraints are shown as horizontal dotted lines in the figures. 

Since x and <p are defined in a model-dependent manner, it is also useful to show the 
corresponding contours on the M z >-M W i plane to compare different G(221) models. We 
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FIG. 1: The 95% confidence contours of the various models on the x-cos(p plane, with Mjj and 
nit either fixed as SM best-fit values (solid) or allowed to be re- fitted (dashed). The area to the 
right of each curve (large x value) is allowed by the global- fit analysis. 
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UU-D NU-D 




FIG. 2: Same as Fig. [TJ but for the models that follow the breaking pattern II. The solid and 
dashed lines are almost indistinguishable. 



translate the constraints on the parameter space of x and <p to constraints on the masses of 
the new heavy gauge bosons, and plot these bounds on the M Z '-M W / plane in Figs. [3] and 
|4j From these plots, we can read off the lower bounds on the masses of the Z' and the W 
in these models, which are presented in Table VII In the UU-D and the NU-D models, the 
masses of the Z' and the W bosons are nearly degenerate, as shown in Fig. |4j In these two 
models, the minimum masses of the Z' and the W consistent with the current experimental 
data are respectively 2.48 TeV and 3.56 TeV. 



TABLE VII: Lower bounds on the masses of the new heavy gauge bosons. The superscripts SM 
or NP indicate whether Mh and fht were set to their SM best-fit values or fitted with rest of new 
physics parameters. Compare with the plots of the Mz'-M\y/ plane in Fig. [3j 





Mf, M) [TeV] 


Af|T [TeV] 


[TeV] 


M*F [TeV] 


LR-D 


1.602 


1.602 


0.269 


0.269 


LP-D 


1.752 


1.742 


0.697 


0.695 


HP-D 


1.674 


1.673 


0.403 


0.403 


FP-D 


1.685 


1.583 


0.673 


0.665 


LR-T 


1.607 


1.607 


0.197 


0.197 


LP-T 


1.753 


1.745 


0.495 


0.493 


HP-T 


1.680 


1.679 


0.289 


0.289 


FP-T 


1.687 


1.587 


0.478 


0.472 


UU-D 


2.479 


2.474 


2.479 


2.474 


NU-D 


3.562 


3.558 


3.562 


3.558 
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FIG. 3: The 95% confidence contours of the various models on the Mz'-Myy/ plane, with Mh and 
Wit either fixed as SM best-fit values (solid) or allowed to be re- fitted (dashed). The dotted lines 
represent lines of constant cos0 and sin</> at fixed values of 0.1 (0.18) for the first (second) breaking 
pattern. Outside the cone surrounded by these regions, one of the gauge couplings becomes non- 
perturbative. 
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FIG. 4: Same as Fig. [3] but for the models that follow the breaking pattern II. The contours 
appear as lines because in these models, Z 1 and W are highly degenerate due to the pattern of 
symmetry breaking. 



B. Key Observables and Their Impacts 

For the models that follow the first pattern of symmetry breaking, we see that models in 
which Higgs triplets break SU(2) 2 x U(l)x lead to smaller bounds on x compared to models 



where Higgs doublets break the symmetry. This is not surprising given Eqs. (41), (46) , 
where we see that, in the triplet models, the corrections to the definitions of the reference 
parameters are suppressed compared to the doublet models. However, the bounds on Mjy 
and Mz> are comparable. 

According to how the contours in parameter space are shaped, the considered (7(221) 
models can be separated into three classes: 

• In the LR, HP, and UU models, large values of eg (sg for UU) are ruled out at small x. 
At small eg (sg), the parameter contours, however, extend to relatively low x values. 

• The contours of the LP and FP models are, by contrast, located at comparatively 
larger x values for small eg. Increasing the values of eg to about 0.8, the contours of 
these models curve to the right (towards higher x). However, increasing eg further 
beyond about 0.8, the FP contours bend towards lower x, while the LP contours 
towards higher x. 

• The parameter contour of the NU model is unique as it is the only curve that bends 
to the left with smaller x value when going up along the vertical axis with increasing 
St value. 

The similarities and differences between the parameter plots can be traced back to certain 
key observables. That is, in the excluded regions of parameter spaces we consistently observe 
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a pattern that several key observables contribute with especially large pulls to x 2 , and it is 
these observables that drive the overall shapes of the curves in Fig. [T} 

TABLE VIII: Overview of the observables that drive the parameter plots. The most and second 
most important observables are respectively marked with the symbols © and ©. In the UU-D and 
NU-D models only one observable significantly contributes to x 2 . 



Model 


Chad 


A FB (b) 


(9l N ) 2 


Qw( 133 Cs) 


Set of other obs. 


LR, HP 


© 


® 








LP, FP 






© 


© 




UU 


© 








© 


NU 






© 




© 



In Table VIII we give an overview of these observables that effectively drive the results in 
Fig. [Tj For each model in the breaking pattern I, we list the two most important observables. 
For the models of the breaking pattern II, we give only one such observable. It is important 
to note that Table VIII only presents qualitative observations that indicate tendencies, and it 
may be the case that some particular points of the parameter spaces have other observables 
that contribute with larger pulls than the ones we indicate here. Nevertheless, the patterns 
we give here are useful in indicating qualitatively which observables are likely to be sensitive 
to new physics contributions from the G(221) models. 

The explicit expressions for the new physics corrections to these key observables, are listed 
in Tables |Lx| and |x| At the 95% C.L., the set of observables (cr had , A FB (b), A LR (e), (g" L N ) 2 '', 
Qvk( 133 Cs)) is respectively measured with a precision at the (0.18, 3.16, 2.80, 0.90, 1.24) per- 
cent level. 

Based on these expressions we can roughly reconstruct the respective shapes of the con- 
tours, and in Fig. [5] we illustrate our argumentation. 

We find that the shapes of the contours for LR, HP and UU models are driven by cr had . 
For the LR and HP models, A FB {b) and A LR (e) also play an important role. Since A FB {b) 
is defined as A FB (b) = ^A LR (e)A FB (b), the large coefficients in xSA b FB /A b FB (See Table IX) 
originate from the smallness of the SM value of A e LR . The combination x5A b FB /A b FB imposes 
about the same constraints as those derived from x5A e LR /A e LR . Since they are strongly 
correlated, we only list the observable A FB (b). In the LP and FP models, with low 
values, Qw ( 133 Cs) is the most important observable, because of the large constant terms 



independent of or (3 in Qw, as shown in Table |IX| The constant term for either the LR-D 
or the HP-D model vanishes because it is proportional to (T| — T R ) of the electron with the 
quantum number assignment T| = T R = — 1/2. With high values, the observable (g" L N ) 2 
determines the shape of the parameter contours. The NU contour is mainly driven by the 
pull of (g v L N ) 2 ■ ( We note that, A FB (b) has a similar effect on the LP and FP contours as 
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TABLE IX: Fractional new physics corrections to the observables a had , Afb(P) and Q 



w 



(133 



Cs) 

relative to the corresponding SM predictions. To obtain the new physics shifts in the triplet versions 
of the 
by \. 



of the LR, LP, HP and FP models, the prefactors of s 2 - need to be multiplied by \, all other terms 





% <5<3"had/ Chad.SM 


x SA FB (b)/A FB:SM (b) 


LR-D 
LP-D 
HP-D 
FP-D 


-1.13 c\ - 0.142 c\ + 0.0432 s 2 ~ a 

<p q> 2p 

+0.346 A - 0.142 c\ + 0.0432 s 2 ~ 

<p <p 2p 

-1.38 c\ - 0.142 c\ + 0.0432 s 2 ~ a 

<p (p 2p 

+0.0985 c\ - 0.142 c\ + 0.0432 s 2 3 

<p <p 2p 


-30.0 c\ + 67.6 c\ - 20.6 s 2 ~ a 

<p <p 2p 

-46.1 c\ + 67.6 4 - 20.6 s 2 ~ 

ip q> 2p 

-30.9 c\ + 67.6 c\ - 20.6 s 2 ~ 

<p <p 2p 

-47.0 c\ + 67.6 4 - 20.6 s 2 ~ 

<p q> 2p 


UU 
NU 


-0.889 s\ - 0.0132 4 
+0.583 s\ - 0.0132 4 


+0.161 4 + 6.29 4 
+14.2 4 + 6.29 4 



LR-D 
LP-D 
HP-D 
FP-D 



UU 
NU 



8Q W ( 133 Cs) /Q W , SM ( 133 Cs) 



-0.855 4 -0.145 s 2 ~ 
+3.35 - 1.95 4 - 0.855 4 - 0.145 s 2 ~ 

<p <p 2p 

-0.855 4 -0.145 s 2 ~ 

4> 2/3 



+2.95 - 1.95 4 



0.855 4 - 0.145 s 2 3 



-0.855 4 
+0.406 + 0.594 4 - 0.855 4 



TABLE X: New physics corrections to the observable {g v L N ) 2 relative to the prediction of the SM. 







LR-D, LP-D, HP-D, FP-D 
LR-T, LP-T, HP-T, FP-T 


0.0875 + 1.91 4 + 0.839 c\ - 2.84 s 2 ~ 

<p <p 2p 

0.0219 + 0.478 4 + 0.210 c\ - 1.42 4s 

4> <p 2/3 


UU-D 
NU-D 


0.839 4 
2.58 - 0.583 4 + 0.839 s\ 



on the LR and HP contours, though subdominant compared to the other observables. The 
same applies to (g v L N ) 2 for the LR and HP models.) 

The starting point of our discussion is that the SM represents the best description of the 
present experimental data, and the G(221) parameters have to be chosen such that they 
minimize the new physics shifts. We first focus on models in the breaking pattern I, the LR, 
LP, HP, and FP models, and note the following points: 
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X 



FIG. 5: Sketches illustrating the influences of some key observables on the parameter bounds for 
the models of the first breaking pattern. The UU-D parameter contour is driven by cr had as well. 
In the NU-D model {g v L N ) 2 is the most important observable. 



The observables A FB (b) and {g" L N ) 2 prefer small c^. In 5(g" L N ) 2 , the c~ and c~ terms 
have the same sign so that they cannot cancel each other. ApB{b) has the largest 
effect on the allowed parameter space in the LR and HP models as the coefficients of 
and c£ are large in magnitude. 

In the case of the LR and HP models the c~ and c~ contributions to 5cr had have the same 
sign and the s\ 1= term is suppressed by a small prefactor. The pull of cr had therefore 
represents the hindrance for the LR and HP models to accommodate large values of 
at smaller x. 

The large impact of Qw ( 133 Cs) on the LP and FP bounds is due to the large constant 
term in 5Qw ( 133 Cs). The only way to make 5Qw ( 133 Cs) small (other than simply 
raising x) is to have large so that the negative contributions from the c~ and en- 
terals can compete with the positive constant term. Consequently, the low-c^ region 
is ruled out in the LP and FP models and the parameter contours start at higher x 
values than in the LR and HP models. 

In the Chad., A FB {b) and (g 1 ^) 2 observables, the and c~ terms always have opposite 
sign. The LP and FP parameter plots suggest that, depending on the exact interplay 
between and c$, the terms may be able to overcome the c~ contributions such 
that the contours are pulled back towards lower x values. Note, however, that the 



expressions given in Table IX cannot explain the branching between the LP and FP 
contours. To account for that effect we certainly would have to extend our discussion 
to a larger set of observables. 

After these comments on the models of the breaking pattern I, it is easy to understand 
the shape of the UU and NU contours. In the UU model all shifts that we present in Tables 
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IX and X favor small values. Especially the fact that the s~ and terms in <5cr had have 
the same sign leads to the exclusion of the high-s^ region. For that reason the UU plot looks 
similar to the plots of the LR and HP models. The contour of the NU-D model is mainly 
influenced by the correction to (g" L N ) 2 . Since 5 (g v L N ) 2 is small if takes some intermediate 
value we observe a bump in the NU-D contour towards lower x values for values around 
0.7. 

An important consequence of identifying these observables is that we may anticipate the 
future impact that the upcoming measurements, with greater precision, may have on the 
global analysis. For example, the Q-weak collaboration [37] and the e2ePV collaboration [38] 
at Jefferson Lab are expected to have ultra-high precision measurements of the weak charge of 
the proton Qw(p) an d the electron Qw( e ) (with a fractional uncertainties of respectively 4% 
and 2.5%). As Qw{ 133 Cs) is a key observable in driving the results for the lepto-phobic (LP) 
and fermio-phobic (FP) models, we would expect that the future measurements of Qw(&) 
and Qw(p) would have a great impact on the global-fit analysis. This is demonstrated 
in Fig. [6j where we perform the global-fit analysis with the expected future uncertainty 
of Qwip) an d Qw( e )- We find that the LP-D contour is drastically different than those 
presented in Fig. [I] As a result of these further constraints from the Qw(p) and Qw{ e ), the 
allowed region in the Mz> — My/ 1 plane shrinks as well. The lower bounds for the W mass 
increase, for instance in LP-D model, from about 0.7 TeV to 1.3 TeV. 



Assumed rtwulls iroio 
upcoming experiments: 



Qw(p) = 0.0715 ± 0.0029 [Qweak coUalx] 
Q w (e) = -0.046900 ± 0.001079 [e2ePV collab.] 
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FIG. 6: The x — cos(0) and Mz'-My/r contours for the LP-D model with an expected uncertainty 
in the Q-weak data and the e2ePV data. These plots should be compared with the corresponding 
plots in Figs. [T] and [3j and demonstrate that, since Qw(e) is a key observable for the LP model, an 
increase of precision in its measurement has a large impact on the global analysis. In particular, 
at low ct, the lower bound on x, is substantially increased. 
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C. Constraints from Triple Gauge Boson Couplings 



Though we do not include the shifts to the triple gauge boson couplings (TGC) in our 
global analysis, they are nonetheless precise measurements at LEP-II that can be used to 
constrain models of new physics. In particular, the ZWW vertex is measured to a precision 
of about 2% and may be used to constrain models of new physics [28]. In this subsection, 
we compute the shift to the ZWW vertex, and use it as a complementary constraint when 
we discuss the results of our global analysis. 

The shift in the ZWW vertex can be parameterized by the Hagiwara's parameterization 

(<?f ) , (93) 

where the SM value of gf is unity. At LEP-II, using a partial waves analysis, the measured 
ZWW vertex is [H] 

gf = 1.001 ±0.027 ±0.013, (94) 

where the uncertainties are respectively the la statistical and systematic uncertainties. 
Adding these uncertainties in quadrature, we have 

gf = 1.001 ± 0.030. (95) 

It is important to note, however, that the measurement of the ZWW vertex is extracted from 
the process e + e~ —> W + W~ using an event shape analysis (including angular distributions) 
assuming SM couplings for all other vertices. Thus, to properly use the experimental results, 
we have to compute the entire e + e~ —> W + W~ amplitude and attribute all the shifts in the 
amplitude to At/f . 

The full amplitude of the tree-level process e + e~ — > W + W~ is the sum of three diagrams: 
s-channel 7-exchange, s-channel Z-exchange, and t-channel i/-exchange. We denote these 
three amplitudes respectively as Ay, Az, and A v . Even though the LEP experiments utilize 
unpolarized e ± beams, it is useful to consider amplitudes with specific helicities for both e ± 
and W and employ the helicity amplitude method jH] to analyze the amplitudes. The 
individual amplitudes are computed in Hagiwara et al. an d here we only note the key 
features of the dependence in the polar angle 9. For a given configuration of incoming and 
outgoing particle helicities, all three amplitudes A^,z,v are proportional to the Wigner's d- 
matrix elements. For the s-channel amplitudes, this is the only dependence in the scattering 
angle 9 (not to be confused with the weak mixing angle) 

A,, z ocdtj A \9). (96) 



For the t-channel, i/-exchange, we have additional 9 dependence from the //-propagator 

C 

1 + /3 2 -2f3cos9 



A^(B- - ■ m ° n ._ Q ) d% AX (9), (97) 
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where B and C depend on the helicity configuration, but are independent of 9. Thus, for 
a fixed helicity configuration of e + e~ — > W + W~ , we can perform a partial wave analysis to 
project out the df^j AA (#) component of the amplitudes 



A 



A tA(T,AA 



(9)d(cos9) 



(98) 



The uncertainties in gf (which we denote as Ag l z ) can be used to constrain the (2(221) 
models by first identifying 

^(221) + ^(221) + jOQM) = ^SM + ^ + Ag Z^ 



SM 



(99) 



and then express each amplitude in the G(221) model in terms of the reference and fit 
parameters 



7G(221) 



1 + -5 A 

x 



(100) 



where the fractional shifts in the amplitudes 8A lt z,p are functions of (f> and /3. For a fixed 
helicity configuration of e + e~ — > W + W~ , we can compute A~, : z,u for both the SM and 
G(221) models and obtain an excluded region on the x-cos0 plane. The regions that are 
allowed by all the helicity combinations are shown in Fig. [7j We note that, generally, the 
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FIG. 7: The regions on the x-cos<p plane excluded by experimental measurements of Agf (in 
grey) for models in the breaking pattern I, for the cases sin 2 2/3 = and 1. 



bounds given by Agf are more relaxing than those obtained from the global-fit analysis 
earlier. In fact, for small values of /3, the regions excluded by Agf are already excluded by 
the global-fit analysis presented in the previous subsection. 
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VI. CONCLUSIONS 



In this work we analyze the constraints on the masses of the heavy gauge bosons of the 
(7(221) (including the left-right (LR), leptophobic (LP), hadrophobic (HP) and fermiophobic 
(FP) as well as the ununified (UU) and non-universal (NU)) models in a unified view based 
on the classification of the (7(221) models in terms of the patterns of symmetry breaking and 
the gauge couplings of fermions. Adapting the framework of effective field theory, we give 
the effective Lagrangians at the electroweak scale and low energy scale , applicable to any 
(7(221) model, and perform a global-fit analysis about a set of 37 electroweak observables, 
including Z pole data, the mass and the width of the W ± boson, and various low-energy 
observables. The experimental precision with which these observables have been measured 
allows us to put strong bounds on the parameter space of the (7(221) models and to constrain 
the masses of the Z' and W !± bosons. At the same time, we show that the bounds from 
the triple gauge boson couplings on the parameters do not affect the result of the global-fit 
analysis. We present our key results in terms of 95% C.L. contours of the allowed regions 
both on the x-coscj) plane, as well as on the Mz>-Mw> plane, from which we can readily 



give the lower bounds on the masses of the W and Z', which are presented in Table VII 



which can be used as a guide for future collider search. We show that, in the first breaking 
pattern, although the mass of Z' is about 1.7 TeV, the mass of the W in some models can be 
relatively light (of a few hundreds of GeV), particularly in the left-right (LR), hadrophobic 
(HP) models. In the case of the second breaking pattern, due the near-degeneracy between 
the masses of the Z' and the W f± , the W is necessarily heavy. 

In addition to the constraints on the parameters and bounds on the extra gauge boson 
masses, we also find associations between certain key observables and the (7(221) models 
discussed in this paper. As these observables are responsible for 'driving' the shape of 
the 95% contour plots, future measurements on these particular observables would have a 
tremendous impact on our results. We demonstrate such an example in Fig. |6j showing 
that an anticipated precision on the measurement of Qw( e ) could largely increase the lower 
bound on the W mass from the current value of 0.7 TeV to 1.3 TeV in the LP-D model. 

In this work, we focus on the interactions of the heavy W and Z' bosons to fermions. 
To extend our results to include flavor-dependent observables, such as the branching ratio 
Br (b — > 57) and the anomalous magnetic moment of the muon, requires a detailed specifi- 
cation of the flavor sectors of the (7(221) models. Though it is difficult to enumerate the 
many models in the literature, in the advent of the Large Hadron Collider (LHC), it would 
be useful to extend to the flavor sector the insights provided in this work. Moreover, the 
direct search of the W and Z' bosons at the Fermilab Tevatron could further constrain the 
(7(221) model parameter space. The potential of the Tevatron and LHC to observe the W 
and Z' bosons will be presented in a separate work. 
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Appendix A: Effective Lagrangians 



In this appendix we show how to obtain the coulping coefficients in the effective La- 
grangian. There are two effective Lagrangians in our framework. The first effective La- 
grangian is SM-like, and is applicable at the electroweak scale after integrating out Z' and 
W, we parameterize this Lagrangian as 



£™ = tf ZJJ Zpf<fPJ + g R ZJf ZJYPnf + [g L wTif W^f 1 YPj2 + g^Wjf^Pnfi + h.c. 

- B ((7i/ 2 ) l ,h > (7 3 /4) L , R ) ^ (7i/2) LR (7 3 / 4 ) L , R • (Af ) 

Except for the non- universal (NU) model, the couplings g^ f , g L ^ f and 

B f(/i/2) LR i (/3A) LR J are an flavor-universal. The second effective Lagrangian is 
Fermi's theory of four-fermion interactions. It is obtained upon further integrating out 
Z and W bosons, and is applicable below the electroweak scale. We parameterize this 
Lagrangian as 

= -c ((7x/ 2 ) LiR , (7 3 /4) L , R ) ^ (7i/ 2 ) L>R (7 3 /4) L , R • (A2) 

The above coefficient functions 9gjp $wTh an< ^ ^ ((W 2 )lh ' (^3^ 4 )lr) can De obtained 
by the following steps: 

• write down the specific form of the effective Lagrangians at the electroweak scale and 
low energy scale by plugging the formulae in Tables IV and [V] into the Eqs. (26) and 

©; 

• extract the coulpings ^wjh an< ^ ^ ^ com P ar i n & the parametrized Lagrangians 
in Eqs. (Af ) and (A2) with the specific form of the effective Lagrangians; 

• replace all the model parameters e, v, and 6 in the expressions of 5^7/ ' 9wjif 2 > ano - 
C by reference parameters a, Gf and Mz, using the relations in Eqs. (46), (45) and 

For future reference, we list below the coefficient functions q h '^ r , Q^J^r, and 

C n/i/2) LK j (/ 3 /4) LH j in terms of the model parameters. We also give some examples 
of the final form of the coefficients in terms of the fit parameters. 



1. The LR-D, LP-D, HP-D, and FP-D Models 

For the four models that follow the first breaking pattern with a doublet (LR-D, LP- 
D, HP-D, and FP-D models), the difference in the coefficient functions is originated from 



the quantum numbers of the fermions. In Table XL we give the quantum numbers of the 
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TABLE XI: The charge assignments of the SM fermions for the first breaking patter. T^ L and 
are respectively the third component of the isospin for the SU{2)\ and SU (2)2 gauge groups (which 
are conventionally called SU (2) l and SU (2)# in left-right models). These charge assignments apply 
to all three generations. 
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fermions, and present the coefficients of the effective Lagrangian in terms of these quantum 
numbers. 

Performing the above procedure, we obtain 



g whh 



S 20 



(for fi t 2 as quarks in LR and LP, and leptons in LR and HP) 



V2s § * 

0, (for f lj2 as quarks in HP and FP, and leptons in LP and FP) 

(A4) 



9z Jf = — \ T l- 4Q + -4rC3 - Q) I > ( A5 ) 



c 2 7 



9 R ZJ f = ^rl- s lQ+~4(T R -4Q)), (A6) 



SnCn \ X 
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and the neutral-current four-fermion coupling coefficients 

„4 2 



\C ( (7/) L , (/'/') J = (1 + | - ^)(Ti - s|Q)(Tf - sfQ') 

+ [(T 3 L - sjQ)(T' L 3 - Q>) + (Tf - Q)(Tf - s?Q')] 



+ ^(T L 3 -g)(Tf -g'), (A7) 

3/ 



^((7/) R ,(77') R )Hi + |-|)(-^)(-^) 

'(T 3 R - slQ)(-sjQ') + (-4Q)(T^ 3 - slQ') 



c 2 7 



+ ^ 

X 

+ 1 ~( T R- s lQ)( T R- s lQ')> (A8) 

„4 „2 



((7/) L , (//')„) = (i + 1 - - 



+ f|| (r 3 _ Q )( _^) + | (T | _ s 2Q )(r « _ S 2Q') 

+ |(Tl-Q)(n 3 -^Q'), (A9) 

where T| R and Q are the isospin charge and electric charge of the fermion /, and T'^ R and 
Q' are the isospin charge and electric charge for the fermion /', respectively. To obtain 
^ ((^)r' (/'/') )' we neec ^ ^° exc hange Q with Q' and T| with in the coefficient 

function C ((//) L , (//') )■ F° r the charged-current four-fermion coupling coefficients, 
we only list the results for LR-D models as follows: 



c((7i/2) L >(7 s /i)J = i, (aio) 
^((7 1 / 2 ) L ,(7 / 3 /l) R ) =^((7 1 / 2 ) R ,(7 / 3 /l) L ) = 8 f, (ah) 

c((7i/0h.(7.^J4 (A12) 

The final form of fl^y/ an d C can be obtained by replacing the model parameters by the 
reference parameters a, Gf and Mz- Below we only list the results for C ((ww) L , (ee) L ) in 
the LR-D model: 

9z* = 9zi: + Wz* = -0-197 + i (-0.348 + 0.348^- + 1.07s? - 0.718sj) , (A13) 

C ((««) L , (ee)J = C SM ((mw) l , (ee)J + ((««) L , (ee)J 

= -0.183 + i (-0.534 + 0.534s^ + 1.50s? - 1.13sj) . (A14) 
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2. The LR-T, LP-T, HP-T, and FP-T Models 



The coefficients of the effective operators in the LR-T, LP-T, HP-T, and FP-T models 
take a similar form as those presented above for the LR-D, LP-D, HP-D, and FP-D models 
respectively, with the following replacements after applying the identity = 1 — s~: 

*i -> \ s % ( A15 ) 

4 -> H (A16) 
s h ^ ( A17 ) 

and terms that are suppressed by 1/x, but without the factors listed above, are further 
divided by a factor of 4. For example, for the LR-T model, we have 

g \_ e = -0.197 + - (-0.087 + 0.174s^ + 0.268s| - 0.180sj) . (A18) 
3. The UU and NU Models 

For the un-unified and non-universal models, we classify the fermions as 

f = leptons 

UU : < , (A19) 

I F = quarks, 

. f = fermions of the third generation, 
NU : { (A20) 
F = fermions of the first two generation. 



We denote the coefficients of the effective Lagrangians with the notations in Eqs. (Al) and 



(A2). Compared to the first breaking pattern presented earlier, there are considerably less 
coefficients because there are no right-handed charged currents in the NU and UU models. 
Similar to the LR-D model, we obtain the electroweak couplings: 

*™ = A i 1 - fj ' (A2i) 
= ^ i 1 + V) • (A22) 

^,__L|-i*_^ 4 _ :/-]. (A2i) 



sift = 3 zff = 7~r (-4Q) . ( A25 > 
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and the neutral-current four-fermion coupling coefficients 



4 



\C ((//) L , (/7')J = (1 + |)(T 3 - 4Q)(T' 3 - s]Q>) 



4 

— ^ ((T 3 — sjQ)T /3 + T 3 (T /3 - s]Q') - T 3 T /3 ) , (A26) 

Ob 

\c ((7/) R , (77')J = (i + f )(-4^)(-*n, (A27) 

4 4 

^ ((7/) L , (77') J = (1 + |)(T 3 - s}Q)(-s}Q') - S jT\-slQ'l (A28) 
\C ((FF) l , (f'f')J = (1 + |)(T 3 - 4Q)(T' 3 - ajtf) 

+ ^ [T 3 (T' 3 - s?g') + (T 3 - s?Q)T ,3 l + 5t 3 T' 3 , (A29) 

1 „ / A=/ „A \ b ,h. 



C ((FF) r , (F'F')J = (f + f)(-s}Q)(-s}Q% (A30) 

|c ((^) L - (f'f') J = (f + |)(t 3 - 4g)(-4g') + ^T 3 (- S fg'), (A3i) 



and 

1 /^.AN 



2 c ((ff) l , (//') J = (i + 1 )(t 3 - 4g)(T' 3 - s fg') 



^4 5^ 5^ 

X y x y X 



/-=/ „A \ ^ 



C ((FF) r , (//')J = (1 + -§)(s}Q)(-s}Q'), (A33) 
|C ((FF) l , (7?) J = (1 + |)(T 3 - S JQ)(- S JQ') + ^T 3 (- S ?g'), (A34) 

4 4 

\C ((FF) L , (77') J = (1 + |)(T 3 - S jQ)(- a jQ0 - |(-^) T ' 2 - (A35) 

We can get C ((//) , (f'F')) by exchanging Q Q' and T 3 T' 3 in the expresion 

C ( (i 71 -^) , (/ /') ) • For the charged-current four-fermion coupling coefficients, we only list 
the results for UU-D: 

C((7i/ 2 ) L ,(7 3 /4) L ) = 1, (A36) 
C((^2) L ,(7 3 /4)J =C((J 1 h) h ,(F 3 F,)J =1, (A37) 

C((F 1 F 2 ) l ,(F 3 F 4 ) l ) = 1+ 1 -, (A38) 
For the final forms in terms of fit parameters, we only list the results for g\ Ze and 
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C {{uu) L , (ee) L ) in the UU-D model: 

Qzee = tile + fee = "0-197 + 0.0227^, (A39) 
C ((«tt) L , (ee)J = C SM ((««) t , (ee)J + 5C {{uu\ , (ee)J 

= -0.183 + i (o.234s^ - 0.424sj) . (A40) 
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